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Q_[ ABSTRACT: We argue that the strongly coupled quark-gluon plasma can be considered 

Jjn ; as a chiral superfluid. The "normal" component of the fluid is the thermalized mat- 

ter in common sense, while the "superfluid" part consists of long wavelength (chiral) 
fermionic states moving independently. We use several nonperturbative techniques to 
demonstrate that. First, we analyze the fermionic spectrum in the deconfinement phase 
(T c < T < 2 T c ) using lattice (overlap) fermions and observe a gap between near- zero 
modes and the bulk of the spectrum. Second, we use the bosonization procedure with a 
finite cut-off and obtain a dynamical axion-like field out of the chiral fermionic modes. 
Third, we use relativistic hydrodynamics for macroscopic description of the effective 
theory obtained after the bosonization. Finally, solving the hydrodynamic equations in 
gradient expansion, we find that in the presence of external electromagnetic fields the 
motion of the "superfluid" component gives rise to the chiral magnetic, chiral electric 
and dipole wave effects. Latter two effects are specific for a two-component fluid, which 
provides us with crucial experimental tests of the model. 
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1 Introduction 

Chiral properties of the strongly coupled quark-gluon plasma (sQGP) have attracted 
much attention in light of recent measurements performed within the heavy-ion pro- 
grams at RHIC and LHC. The analysis of charge-dependent azimuthal correlations by 
STAR [1, 2], PHENIX [3] and ALICE [4] collaborations suggests a possible local V- 
violation in strong interactions, which manifests itself as an asymmetry in the charged 
particle production with respect to the reaction plane. One of the most popular theo- 
retical approaches to study the observed phenomena is to seek for new electric currents 
of a specified direction in the QGP phase (see e.g. [5-7]). Since the physics of sQGP 
is essentially nonperturbative, there is a lack of models based on the first-principle 
calculations, i.e. starting from the QCD Lagrangian. In our paper we establish such 
a model for sQGP in the range of temperatures T c < T < 2 T c , i.e. slightly above the 
deconfinement transition. We choose the range mainly because of two reasons: first, it 
is estimated to be typical for sQGP at RHIC [8]; and second, we can use hydrodynamic 
models to describe the system [9-15]. 
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Lattice calculations [16-18] demonstrate that the Dirac spectrum for massless 
quarks at these temperatures contains a peak near zero virtuality separated by a gap 
from the bulk of the spectrum (Fig. 1(b) and Fig. 2). Therefore, it is natural to in- 
troduce a two-component fluid model for the sQGP: one component, carrying chiral 
properties of the fluid, and the rest, corresponding to the bulk of the spectrum. In 
our case it turns out that the fluid is described by a system of equations similar to 
ones of a relativistic superfluid [19, 20]. This fact, of course, does not directly lead 
to the conventional superfluidity, since the "normal" and "superfluid" components in 
our case carry different U(l) charges. We also do not assume any hidden symmetry 
spontaneously broken in the system. Instead, we treat the term "superfluidity" in a 
phenomenological sense traced back to the Landau's formulation [21], i.e. a combina- 
tion of two independent (curl-free and "normal") motions of the fluid separated by an 
energy gap. 

Before proceeding with a formal derivation let us provide a few hints supporting 
the existence of the "superfluid" component. In Refs. [22-27] it has been discovered 
and confirmed by [28, 29] that topological charge density forms long-range coherent 
global structures around a locally one-dimensional network of strong fields (so-called 
"skeleton"). Though the simulations were performed for low temperatures, the skele- 
tons could very well survive slightly above the deconfinement transition, since there 
is a corresponding long-range order of gluonic fields [30] in the case of nonvanishing 
topological susceptibility. The extended character of the skeleton can be interpreted 
as a long-distance propagation of (quasi-) particles, carrying finite chirality and thus 
forming the chiral "superfluid" component. 

In addition, the lattice data [22-29, 31, 32] suggest that the topological charge 
density itself (for uncooled configurations) is localized on low-dimensional defects with 
fractal dimension between 2 and 3, i.e. presumably on (percolating) central vortices (see 
also [33, 34] for a similar result based on the scalar density distribution and [35, 36] for 
a localization on vortex intersections). One can demonstrate [34, 36-38] that removal 
of the central vortices eliminates all zero and chiral near- zero modes from the fermionic 
spectrum. The converse statement is also true: small eigenvalues can be generated at 
separate center vortex configurations [37]. It seems natural to consider the vortices as 
spatially one-dimensional "guides" for light fermions. Since massless fermions propagate 
in (1+1) dimensions with the same speed (speed of light), we can consider their bilinear 
combinations as bosonic excitations, even without assuming their interaction with each 
other. A similar effect of binding for light fermions along ID defects, without making 
use of the Goldstone theorem, was considered in [39, 40] for non-abelian fields and in 
[41, 42] for a strong external magnetic field. 

It is also known that monopole trajectories and central vortices populating QCD 
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vacuum and leading to the confinement at low temperatures [43] become Euclidean- 
time oriented (static) at high temperatures [44]. Static nature of the vortices makes it 
possible to continue them to the Minkowski space-time by considering their intersections 
with planes of constant Euclidean time. These intersections are in general percolating 
strings, which, according to the polymer representation of the field theory [45, 46], can 
be interpreted as a 3D scalar field with non-vanishing vacuum expectation value. If this 
field is complex, then its phase can describe a new Goldstone mode in the deconfinement 
forming the superfluid component of sQGP [47, 48]. Taking into account the vortex 
picture of the topological charge generation [49, 50] it seems interesting to apply the 
arguments of [47, 48] to the chiral superfluidity (we do not implement it here). 

Other arguments in favor of the superfluidity were presented in recent analysis in 
the framework of stringy models [51-53], where a new formulation of the superfluidity 
with vanishing chemical potential has been suggested. 

In what follows we do not specify the way quarks are bound to each other and form 
a "superfluid" component, leaving the question open for the further studies (possible 
examples are given above and in the main text). Regardless of the nature of the 
"superfluid" component, we obtain the same universal phenomenological predictions, 
namely the chiral magnetic, chiral electric and dipole wave effects. 



2 Derivation of the effective Lagrangian 

It is known that many of the essential properties of the QCD vacuum (such as the value 
of the chiral condensate, electric and magnetic characteristics, etc.) can be determined 
from the IR part of the fermionic spectrum [28, 54-58]. This makes it possible to 
introduce a finite cut-off A for the fermionic spectrum, without affecting the values 
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of the observables (indeed, our phenomenological results do not depend on A, see 
Section 3.2). In this section we perform the bosonization procedure with a finite cut-off 
[59] for the SU(N C ) x U em (l) theory and derive an effective Lagrangian. This procedure 
leads to appearance of a dynamical axion-like field 6(x), which we identify with the 
propagating chirality (i.e. local difference between numbers of left- and right-handed 
fermions) in the system. 



2.1 The functional integral 

The gauge fields of the theory are represented by 

A, = A°T° + gG;T a = AlT\ A 5fl = A^T° , (2.1) 

where |T a a = 1, (N% — 1) j are the SU(N C ) color matrices normalized by Tr(T a T b ) = 

5 ab /2 and T° = 1. Here A ^ and A^ play a role of abelian fields, while are gluonic 
fields. The axial-vector field A^^ is an auxiliary external field and will be turned off at 
the end of the procedure. The Euclidean functional integral for Dirac fermions (Nf = 1) 
in external vector A^(x) and axial- vector A^^x) fields is given by 

Z(A,A 5 ) = J ' DijiDif) expj- J d A x^{lJ) - im)^J = det(# - im) , (2.2) 

ip = + A + 

The eigenvectors of Ip are defined by 1 

Jj)^ n = A n ^„ . (2.3) 



They form a complete orthonormal basis {"0n|^ G N} in the space of square-integrable 
spinors. According to the general prescription [59-61], we expand ip and jf> into this 
basis and cut the summation after first A" basis vectors: 2 

N N 

ip(x) ->■ anipn(x) , i>{x) -> Y i>H x ) h n • (2.4) 

n=l n=l 

1 The Dirac operator Ip should be Hermitian, therefore formally the axial field has to be rotated 
A$ —> iA$ from the beginning and then rotated back A$ — > —iA$ in the final result. We stick to 
conventions of [59] and choose Hermitian gamma-matrices 7^ = 7 and 75 = —70717273, while in 
Fujikawa's earlier works [60, 61] they are anti-Hermitian. 

2 We do not consider zero-modes, otherwise the functional integral (2.2) will vanish in the chiral 
limit. We are not interested in the global topological properties of the system. 
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Then the functional integral (2.2) can be rewritten as 

Zn = / ] \ db n da n exp I - ^ b n {ip n \lj) - im\ij) k )a k > = det($ 

^ n=l L n,k=l ) 

N 

{$ - im) N = 1-P N + P N ($ - tm)P N , Pn = J2 MM ■ 



im) N , (2.5) 



n=l 



From the gauge invariance of {JJ) — im)^ it follows that the projector P/y commutes 
with the Dirac operator P\r] = 0, which will be used in the further calculations. 

For the regularization we need to introduce a mass parameter A, to be discussed 
below. It will be defined via |Ajv| < A < | A I - The projector P/y can then be replaced 
by 



Pa = 6 1 



A 2 



+00 



27ri((-ie] 



exp 



iC(l-^ 2 /A^ 



;2.6) 



where 6(x) is the Heaviside step function. We also replace all the indices N by A. 
2.2 Vector currents conservation 

The functional integral (2.2) is invariant under the gauge transformation 

-> Af, + d^a + [Ay., a] , 
A§y — > A 5/1 , 



ip — > (1 — a)ip 



(2.7) 



where a = a a T a . Therefore, 



OA^X) 



X8) 



After integration by parts we get 



(2.9) 



which means that the vector currents 

i"*(x) = 1 6Za{A) 
Ja K ' ~ Z \ 5A*(x) 







im A 



| a;) 



(2.10) 
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are conserved. To derive the last expression one can use the equivalences exp Tr ln(-) 
det(-) and 



1-Pa + 



Pa 



im 



(2.11) 



(0 - im) A 

as well as the general identity for projectors 

pp'p = P(P 2 )'P = IPP'P = , (2.12) 
where the prime denotes a derivative of P with respect to an arbitrary parameter. 
2.3 Anomaly for the axial current 

Let us perform the chiral transformation of the functional integral 



A^ — y A^ , 

K K + <V . 

-> (1 + 7se)^ , 
ip -^^(l + 7 5 e), 



(2.13) 



with e = e°T°. Then in analogy with the previous section we get 



Z A (A 5 ) = Z A (A 5 + de, e) ~ Z A (A 5 ) + J d A x 



SA 5 °Jx) 



■ d u e(x) + / d A x 



SZ A (A 5 ,e) 



8e(x) 



e(x 



(2.14) 



The second term here is due to a nontrivial transformation of the integration measure 
of the functional integral, and was absent in the case of vector current. The projector 
Pa arid combined operator ij) — im also transform under the chiral rotations: 



P A ^(l + 75 e)P A (l + 75 e), 
^M = { 75 ,P A (x)}5(x-x') 

If we introduce the axial current as 



($ - im) -)• (1 + 7se)($ - im)(l + j 5 e) 
S(lp — im)(x) 



then from (2.14) it follows that 



5/i 



5e(x') 

1 5Z A (A 5 ) 
Za5A 5 1(x) 

1 SZ A (A 5 ) 
Z A Se(x) 



—2im^<~,5(x — x 1 ) 



(2.15) 
(2.16) 

(2.17) 



(2.18) 
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The latter is equal to 
1 5Z A (A 5 ) 



Z A 8e(x) 



Tr 



10 



im) A 



■ 2^($-zm)P A + P A 



5(^-1 



im) 



Se 



Se 



Se 



-Pa 



2imTr 7 5 T u (:c 



P 



A 



zm) A 



x) +2Tr(7 5 (x|P A k)) 



(2.19) 



where we used (2.11, 2.12, 2.15, 2.16) and the commutation properties of JJ). By means 
of (2.6), the second matrix element in this formula can be expressed as 



Tr( 7 5 (x|P A |x)) = Tr { 7 5 



d 4 k 



+oo 



d( 



(2vr) 4 J 2m(C-ie) 



e -ikx e iC(l-Ip 2 /A 2 ) e ikx 



(2.20) 



d 4 k 



(2tt) 4 J 2m(C-ie) 



cxp 



(2.21) 



d 4 k 



o 



(2tt) 4 J Am(C,-ie 
1 

A2 



75 (# + 4A 2 (PP M ) 2 + A 2 ^ 2 ) 



(2.22) 



where k = fc/A. Using the algebra of the gamma-matrices, we write 

IP 2 = -D 2 V - ^ 7 [ V ] A- + 4£ - (^^5 M )7 5 - 7 [ V ] (£>VmA, - A^^Dyy) 7 5 



(2.23) 



where Py = —i(d + A), and finally obtain 



e^Tr [ A^A XK + - A^A^ 



4vr 2 



Tr du&'&'A 



{{d*As, A Su ] , A^} + i {8M S(1 , A 5 "^ 5 4 + l^ff-A^s^ 



where the chirality p A is defined by 



(2.24) 



p A = zTr ( 7 5 r°(x 



Pa 



«m A 



| a;) 



(2.25) 
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and the field strengths are 



A^ u — A v — dyAp + [Ap, A v \ , 

Mpv = dpA 5u - d v A b p + [An, A 5v ] - [A v , A 5 J = d^A 5u - d v A hii . 

The first line in (2.24) is the usual chiral anomaly with a mass- dependent term and 
the topological term. The rest of (2.24) can be subtracted away by adding to the 
Lagrangian the following gauge-invariant term: 

A£ = ^Tr (3A 2 (A 5fl A 5 ») + 1 {d.A^f - OW) 2 ) - (2.26) 

2.4 Axionic Lagrangian 

The external axial-vector field A^^ should be now switched off, which is equivalent to 
considering the pure gauge A^ = 0^9, since one can always generate a pseudoscalar 
field 9 by the local chiral rotation Ji — > ^4 + rf$9. In other words, the quadratic term 
in (2.26) is required by having the chiral transformation consistent with the correct 
form of the chiral anomaly. In the chiral limit m — > the total effective Euclidean 
Lagrangian is given then by 3 

4 4) = AC + \g^G% + \f^F^ - j^Al - f»A 5 l 

= AC + hcr^G% + \f^F, v - f»Al - gj^Gl - 3 fd,9 

= \G a ^G; u + \f^F, v - f»Al - gf»G« 
A 2 N a 2 N 

+ ^ eu2e ~ V™ 9 * 6 ? ' (2 - 27) 

where we used the anomaly expression (2.24) and integrated by parts. Kinetic terms 
for the electromagnetic and gluonic fields could be introduced already in (2.2) and do 
not affect the derivation. As we see, the kinetic term for 9 even being absent in (2.2) is 
generated dynamically, see also [62] for a similar discussion. If we drop quartic terms 
and replace conventions A ^ — > Ap, j° M — > j^, then the effective Lagrangian is reduced 
to 



4 w ' 4 

A 2 A^ C QU/1Q „ , g 2 , K 

?7T 



3 Currents j afi can be sourced by fermions from the scales above A. 
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This Lagrangian describes a generalization of the axion electrodynamics [63], where 

the new terms are due to gluonic fields. It seems interesting that the axionic field 

6(x,t) appears within QCD coupled to QED, without any further assumptions as e.g. 

the Peccei-Quinn mechanism [64]. Similarly to the 'true' axion, 8(x,t) is a propagating 

2A , — 

dynamical field. However, the value of the decay constant / = — y N c turns out to be 
of order of scales appearing in QCD (see below), while in cosmological scenarios this 
value is usually around 10 9 — 10 12 GeV. The formal similarity of (2.28) to the axion 
Lagrangian allows us to derive an explicit expression for the mass of 8(x,t) [65-67], 

m 2 J 2 = X (T) , (2.29) 

where x = h m v^oo is the topological susceptibility related to fluctuations of the 
topological charge Q. Lattice simulations demonstrate that \ goes (almost) to zero 
at temperatures above the deconfinement transition [16, 17, 68-70], this behaviour is 
also confirmed within the interacting instanton liquid model [67]. Meson masses in the 
deconfinement interpolate between their values at T < T c and approximately twice the 
lowest quark Matsubara frequency (i.e. 2ttT). These two facts allow us to consider the 
axion-like field 9(x, t) as a nearly-massless field, an essential requirement for a superfluid 
mode. 

Instead of considering 9 as a real particle we rather tend to interpret it as a collective 
excitation (quasi-particle) of the medium in the nonperturbative regime of QCD (e.g. 
combinations of chiral quarks). As will be shown in next sections, the excitation carries 
chirality and can be considered as a 4D generalization of the Chiral Magnetic Wave 
[41]. These excitations can also be exactly massless (compare with phonons or sound 
waves) and, at the same time, not necessarily consisting of massless quantum particles, 
appearing as Goldstone bosons of some broken (e.g. Peccei-Quinn-like) symmetry 4 . A 
straightforward derivation of this collective solution made out of quarks and gluons 
is not worked out at the moment. However, there are some evidences that it could 
exist, see e.g. comments on binary bound states in sQGP [71, 72], lattice results on a 
screened attractive force in the color-singlet channel [73] and ideas mentioned already 
in the Introduction. 

2.5 Interpretation of A 

The scale A can be studied by considering Nf = N c = 1 in the limit of a constant 
background /x = const, fi 5 = const and negligible anomaly (i.e. slow varying 9 = i/^is). 

4 Even if we assume that U(1)a symmetry is broken spontaneously, as suggested in [74]. then the 
bosonization procedure does not describe t/, since it is a spin-0 particle, which can not carry chirality 
as 8 does. 
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In this case, keeping also the quartic terms in (2.27) for generality, we get 

•(4) , / A \ 2 



P5 



hm £ 



Pb 



3tt 



2^5 



(2.30) 



In other words, the value of A can be read off from the dependence of chirality ps on 
the chiral chemical potential /15. We consider here three existing examples one can find 
in the literature. 

(1) At high temperatures, neglecting effects of gluons and assuming equilibrium, one 
can define the thermodynamic grand potential as [6] 



8=± 



d 3 p 



UJ 



p,s 



T^log(l + e- 



T ) 



(2.31) 



where w 2 s = (p + sps) 2 + m 2 and one also assumes an approximate conservation 
of the axial charge. Differentiating the grand potential with respect to /i 5 and 
taking the massless limit one obtains [6] 



P5 



2 , P 



71- 



P5 



3tt 



2 /4 ■ 



(2.32) 



Comparing (2.32) with (2.30) we conclude 

^2 



A 



7T 



T 2 + 



7T- 



(free quarks) 



(2.33) 



In [6, 7] this scale is compared with the inverse radius of a typical sphaleron at 
given temperature. Notice, that in the high temperature limit T>/j and T ^> p 5 
we get simplifications A ex T and ps oc A 2 ps oc T 2 \i§. 

(2) In case of a strong external magnetic field (eB > /i 2 /2) one can construct the 
grand potential for fermions on the lowest Landau level [6] 



n 



eB 
4?? 



d 3 P\ 



e T 



± 



(2.34) 



where w 2 = (p|| + P5) 2 + ^ 2 and p\\ denotes a component of momentum parallel 

eB 

to the magnetic field. This gives us ps = — -U5 and hence 

2ir 2 



A = 2\/ eB (free quarks and strong B) 



(2.35) 



- 10 - 



Upon the redefinition 9 —> 9 the kinetic term for 9 in the effective 

y/2N c eB 

Lagrangian (2.27) takes a canonical form —(d^) 2 , while the quartic terms are 

suppressed by factors 1/B and 1/B 2 , respectively. Therefore, the bosonization 
procedure becomes exact in the limit B — > oo as in [41]. 

(3) To include effects of gluons one needs to perform a lattice calculation with finite 
/is, which has been done in [75, 76]. Slope of the curve ps = ^5(^5) obtained 
in the paper 5 is approximately one (in lattice units), which being translated to 
physical units and compared with (2.30) gives us 

A ~ 3 GeV (dynamical lattice fermions, Nf = 2, N c — 3) . (2.36) 

Appearance of this scale (much larger than Aqcd) is not surprising, see e.g. 
[77, 78]. 

It is worth to mention that all three predictions even if being affected by either rought 
initial assumptions or lattice artifacts, still provide a finite and reasonable value of A. 

2.6 Fermionic spectrum and chirality 

The Dirac spectrum for massless fermions is schematically shown in Fig. I. 6 Lattice 
studies can be found in [16-18, 34, 79]. Below the critical temperature T c the spectrum 
is a continuously growing function in A. 

In contrast, at T c < T < 2 T c the spectrum consists of three parts: exact zero modes 
(and near-zero modes), followed by a gap for the low-lying modes and a continuous 
spectrum starting from a finite A#. Presence of the near- zero peak can be interpreted 
as a manifestation of a (small) remaining chiral condensate [16, 18], which is, howewer, 
not yet rigorously proven. Both the chiral condensate and topological susceptibility 
X(T) are defined from the near-zero modes, since the exact zero modes do not survive 
in the thermodynamic limit [55, 69, 70]. It is important, that x{T) is small enough (to 
keep me small), but still not zero (otherwise 9 itself does not exist). This forces us to 
choose the window of temperatures T c < T < 2T C . At higer temperatures T > 2T C 
the peak disappears completely and all the correlations between quarks supposed to 
be washed out by thermal effects. The gap width seems to be temperature dependent 

There is no a priori introduced UV-cutoff in the paper since the inversion of the Dirac operator 
is done by means of the BiCGstab solver 

6 We show here only the non-negative part of the spectrum. In the chiral limit it is symmetric with 
respect to the reflection A — > — A. 
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Figure 2. Typical fermionic spectrum in the deconfined phase as seen by a SU(3) quenched 
lattice simulation with tadpole-improved Liischer-Weisz action and overlap fermions (/3 = 8.45, 
a = 0.095 fm, V = 16 3 x 6, N f = 1, T c ~ 260MeV). 



and grows with the temperature [17]. The right bound Xb is natural to identify with 
an effective quark mass, since on the corresponding fermionic mode 

X 2 b 4j b = fy B = f^ B = m 2 eff ^ B ■ (2.37) 

A very strong external magnetic field can slightly shift the right bound of the spec- 
trum to the left 7 (see Fig. 2), but for magnetic field strengths occuring in heavy ion 
collisions[7] the principal shape of the spectrum remains the same. 

Let us demonstrate that chirality is determined by the first part of the spectrum. 

Pa = * Tr 7 (x — —-\ x )\=i > : — 

v v ^ ' A 7 0<|A|<A 

. ^aV^a , . ^aVtW'a ^aV^a 

= 1 / : 1" 1 > \ : = -2m > -5— - 

*— ' A — im z — ' —A — vm X 1 + m z 

0<A<A 0<A<A 0<A<A 

A 

/2m 
dXu(X)——^i>xh 5 i>x, (2.38) 
X z + m z 



where u(X) denotes the spectral density for an eigenvalue A. Here, as before, we dropped 
the exact zero-modes. Then, using the identity 

777 

lim n xT— - = ^( A ) ( 2 - 39 ) 

m->0 A + 777 



7 We are grateful to Victor Braguta for making this observation from our data 
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we find the following expression for the chirality 

A 

lim p\ = -2tt I d\ z^AWA^aWa = -tt lim i/(A) ^aV^a • (2.40) 

m->0 J A->0 



From this expression we see that the chirality can be determined exclusively from 
the near-zero fermionic modes. It is suggested and tested on a lattice [16] that this 
part of the spectrum is originated from the zero modes of separate topological defects 
populating the vacuum, because interactions between the original zero modes break 
their degeneracy. So the exclusion of the exact zero modes from our analysis does not 
affect much the main results (for additional arguments see also [57]). 



3 Quark-gluon plasma as a two- component fluid 

We now consider Minkowski version of the effective Lagrangian (2.28) for the quark- 
gluon plasma with one quark flavor 8 , 



mdfi + -f9G a ^G% + -^9F^F^ } (3.1) 



2A , — 

where, again, the decay constant is defined by / = — y N c and the anomaly coefficients 

2 7F 

are given by C 7 = ^ and C g = J^. 

Varying the Lagrangian (3.1) with respect to 9 and vector fields and we 
obtain the following equations of motion: 

d"d,9 = -^F^F, V - ^- 2 G^ a G% , (3.2) 

d^^-f + C^d^F^, (3.3) 
8^ ua + gf abc G b ^ luc = -gf a + C g [d K 9)G Kva , (3.4) 
where f abc are the structure constants of SU(N). Bianchi identities are given by 

d li FK u = t (3.5) 
d^ va + gf abc G\G^ c = . (3.6) 

Interestingly, 9(x,t) obeys a wave equation which is sourced by the U(1)a chiral 
anomaly. This makes it a four- dimensional generalization of the so-called "Chiral Mag- 
netic Wave" which has been recently proposed in [41]. 



*The metric we use has the signature ( — h ++) 
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3.1 Hydrodynamic equations 

The hydrodynamic equations may now be derived from the effective Lagrangian (3.1) 
and the corresponding equations of motion (3.2)-(3.6). Taking the divergence of (3.3) 
and using the Bianchi identity (3.5), we obtain the conservation law for the electro- 
magnetic current, d v j u = 0, provided that the topology of the #-field is non-singular 
with 9 

[d ll ,d v ]9 = 0. (3.7) 
Varying with respect to the axial-vector —d^O, we obtain the axial current 

n = fd^6 . (3.8) 
This current satisfies the anomaly equation 

« = ~^% v - ^G^G% , (3.9) 

as can be seen by substituting the equations of motion (3.2) into the divergence of 
(3.8). Notice that U(1)a is still broken in the deconfinement [69, 70]. 

The total energy-momentum tensor is a sum of both electromagnetic 0^" and 
gluonic ones Q£ v and the stress-energy tensor of the fluid T^ v . The energy- momentum 
tensors of the free electromagnetic and gluonic fields are given by 

0f = F^F\ - ^g^FapF^ , (3.10) 

Their divergencies can be found by means of the equations of motion (3.3, 3.4), 

dpQt? = -F» x 3x - C.F^F^e , (3.12) 
d^&g = -gG uXa f x - C g G uXa G XK a d K 8 . (3.13) 

Substituting this into the conservation law of the total energy-momentum, 

d^T^ + 0f + 0f ) = , (3.14) 

we get 

= F ux Jx + gG^aja + Cl F vX F XK d R B + C g G vXa G a XK d K 6 . (3. 15) 



9 delta-function on the r.h.s. of (3.7) would generate an axionic string [80]. For relevant phenomeno- 
logical consequences see [81, 82]. 
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0m# = —i-F^r ~ -fG" va Gl u , (3.18) 



The first two terms on the right hand side are the standard terms for a work done by 
external fields. The last two terms is a work done by a Witten-like current [83]. Finally, 
keeping that p§ = j® = f 2 d°9, we obtain an expression for the axial chemical potential, 

5C {2) 

p 5 = ^— = -d 9. (3.16) 

When boosted this turns into the Josephson-type equation p^ = —u^d^O, where u yi will 
be chosen later as the velocity of the "normal" component, normalized by the condition 
u^u^ = — 1. We also should assume, that d t 9 is slow varying in time, i.e. 8^9 ~ 
0(p°), so the changes in chiral (axial) charge are small at the scale of QGP lifetime 
[84, 85], otherwise the chiral chemical potential is not well-defined. Hydrodynamic and 
constitutive equations are then of order 0(p 2 ) and 0(p x ), respectively. 
In summary, the system of hydrodynamic equations is given now by 

dfjT^ = F v \j x + C.h^O) + G vXa {gf x + C g G XK a d K 9) , (3.17) 

T 

d l j» = d ft j> M = 0, (3.19) 
u^e + (i 5 = , (3.20) 

where the last one is the Josephson-type equation. Corresponding constitutive relations 
in gradient expansion, satisfying (3.2) and (3.17-3.18), can be represented by 

= (e + P) uV + Pg" v + f 2 d"9d u 9 + , (3.21) 

f = pu » + , (3.22) 

f a = p a u » + v^ a , (3.23) 

jg = fd*e + v£, (3.24) 

where e, P, p, p a are the energy density, pressure, electric charge density and color 
charge density, respectively. Terms r^ u , i/ 4 , v^ a and denote higher-order gradient 
corrections and obey the Landau conditions 

%r^ = 0, tt/ = 0, v"" = 0, «/ = 0. (3.25) 

The stress-energy tensor T^ u consists of two parts, an ordinary fluid component and a 
pseudoscalar "superfluid" component. This modifies the Gibbs relation 



dP = sdT + pdp — f 2 d 



2 ' 



(3.26) 
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where s is the entropy density. Being additionally supported by the results of the 
Section 2.6, we can describe the fluid content as a mixture of two components, originated 
from different parts of the fermionic spectrum, 

(1) Zero and near-zero fermionic modes, which are involved into a potential (curl- free) 
motion of the chirality described by = f 2 3^9 = P5^5 M , where us^ = P5 0^9 is 
the "superfluid" velocity. 

(2) The rest forming an "axially-neutral" component described by an electric j M = pu^ 
and color j a = p a u^ currents and separated from the curl-free component by a 
finite gap. 

There are numerous studies (e.g. [7, 86, 87]) suggesting generation of a finite chirality 
in the processes of heavy ion collisions, so we can assume (ps) = / 2 (/is) 7^ within an 
event. Such kind of initial conditions is not captured by lattice simulations ((/O5} = 
after averaging over gauge configurations [57, 58]), unless the finite (/i 5 ) is introduced 
ad hoc [75, 76]. 

3.2 Phenomenological output, possible tests of the model 

Electric and magnetic fields in the fluid rest frame are defined as 

E» = F^u u , B» = F» u u v = -e^ a/ ViV (3.27) 

One can also rewrite these definitions in the following way 

= e^ a ,]U a BP + u»E v - u v E^ (3.28) 
= e^u a E p + u, L B v - u u B^ , (3.29) 

which we will use later. 

The second term in first brackets (3.17) or equivalently the second term on r.h.s. 
of (3.3) has an interesting phenomenological interpretation as an additional electric 
current, induced by #-field, i.e. assosiated with the "superfluid" component. This 
current is conserved due to the Bianchi identity and vanishes in absence of either 
external electromagnetic fields or 9. Let us split this term in three pieces using (3.29): 

jf = CryF\ K d K 9 = C y d K 9 (u K B x + e KXal3 u a E? - u x B K ) (3.30) 
= -C> 5 5a + a ( e XaK pu a d K 9E? - C jUx (d9 ■ B) . 

The first term in the sum is nothing but the Chiral Magnetic Effect (CME) [6], i.e. gen- 
eration of the electric current along the magnetic field. The second term is analogous to 
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the Chiral Electric Effect (CEE) [88], i.e. generation of the electric current perpendic- 
ular to applied electric field and to both (normal and superfluid) four-velocities. The 
third term is a dynamical realization of the domain wall polarization [89, 90] or simply 
the fact that the "would-be" axions acquire an electric dipole moment in a magnetic 
field [63]. Keeping in mind a wave-like profile of the ^-field we could call this effect 
the "Chiral Dipole Wave" (CDW). First term contains /x 5 which is in our case derived 
within QCD in contrary to other CME models. Last two terms are specific for the two- 
component fluid model and can be considered as a possible experimental test of the 
model. A concrete quantitative prediction and comparison to the experimental data 
are beyond of the scope of this article and will be presented in a future publication. 

The same analysis can be applied to the color current 



which is the second term on the r.h.s. of (3.4). The corresponding effects can be called 
"Color CME", "Color CEE" and "Color CDW", respectively. We are not focusing on 
these effects, since the color currents cannot be observed directly. It is also important 
to mention that both currents (3.30, 3.31) do not depend on the UV cutoff A introduced 
in the preceding sections. 

3.3 Change in entropy and higher order gradient corrections 

In previous sections we kept higher-order corrections r^ u , z/ M , v^ a , vt* undetermined. 
These corrections incorporate possible dissipative effects (in presence of e.g. viscosity 
or electrical resistivity) and can be found from the constraint > on the entropy 
current s M [91]. A priori it is not obvious whether they can interfere with the result 
(3.30) or not. Indeed, within the ordinary hydrodynamics the chiral magnetic effect 
can be found as a part of the term [92-94]. In this section we show that there are 
no higher order corrections to the electric current arising in our case and affecting 
the phenomenological results of Section 3.2. Following [92], we transform the quantity 



■Sa _ 

J\ = 




(3.31) 



9 



(3.32) 



using hydrodynamic and constitutive equations and equate both resulting expressions. 



Constitutive equations 

Using the second law of thermodynamics e + P = Ts + pp and (3.21)-(3.24), we can 
rewrite this quantity as 



I = ujQ^ ((Ts + pp) m m ) u v + Uy {Ts + pp) u^d^u" + Uyd^Pg 

2 Fi y f) Fi Fi^fi _l „ f 2 F)VQ Ft F)VQ i „ Ft h-W 



[IV 



T' 



+ u v fd v 6 d^6 + uj'd^e d v d»0 + u v d l 
+ pd^ (pu" + i/*) + psfd^O + . (3.33) 

Using the normalization condition u^u^ — — 1, the Josephson equation u^d^O + /is = 
and the following identities 

<9X = inv = d u° = 0, (3.34) 

UyvPdpU 1 ' = - {u v u»dX + u'vPdnUv) = u»d^{u u u u ) = 0, (3.35) 

d^puT) = u^p + pd^ = u^d^p = inv = u°d p = (3.36) 

we can simplify the expression for / to 



I = — Tdp {u»s) - {sd^T + pd^p - fd u 6 8^9 - d^P) + u u d,,T 



[IV 



+ pdX + p 5 dX + 2C,p 5 E x B x - -^p 5 G^ a G% . (3.37) 

The sum in the curly brackets is identically zero due to the thermodynamic relation 
(3.26). Also Uyd^T^ = —T^d^Uy because of the Landau frame condition (3.25). This 
leads to the further simplification, 

(3.38) 

Hydrodynamic equations 

Let us rewrite / again using the hydrodynamic equations (3.17-3.18). Then 

/ =F vX u v pu x + u v F vX v x + gG vXa u vP a u x + gu v G vXa v a x 
+ C lUv F vX F XK d K e + C g u u G uXa G a XK d K 6 

+ C,p 5 E x B x - 9jL^ a Gl v 
= - E x u x + 2C,p 5 E x B x - ^Lp 5 G>^ a G« u , (3.39) 

where we used the Josephson equation, the definitions of the electric and magnetic 
fields (3.27) and the corresponding inversed relations (3.29). 



- 18 - 



Combining this result with (3.38) we obtain 

-T [d, (su») - - ^dX) = r^d,u v - E x u x . (3.40) 

Then dividing by — T and adding —v^dj^ — v^dj^ to the both sides we can rewrite 
the result as 

% <y - ^ - y v ") = -f( d ^ u - ^ ( 9 4 - t e ») - v " 9 »y ■ (3 - 4i) 

Comparing this result with one of Son and Surowka [91], we see that the term 
proportional to C 1 is absent. This fact tells us that the divergence (3.41) is well 
defined, i.e. the entropy production is always nonnegative, and in contrary to [91] we 
do not need to add any additional terms to the entropy current. Therefore, there are 
no leading-order corrections to jf coming from the dissipative term z/ M . So, all three 
anomalous effects are present already at the level of the hydrodynamic equations (3.17). 

The expression for the entropy current remains unchanged also because the "su- 
perfluid" component itself has zero entropy. Indeed, considering (3.41) in absence of 
dissipative corrections we obtain 

^(sO = 0, (3.42) 

i.e. only the "normal" component cotributes to the entropy current, which is a common 
property of a real superfluid. This fact would well agree with the long-range coherence 
of the chirality distribution[22-29, 95, 96], but should be studied more carefully, since 
the microscopic nature of 9 is not known precisely. 

4 Conclusion 

In this paper we provided a novel treatment of the sQGP dynamics for the tempera- 
tures typical for recent and ongoing heavy-ion experiments. The main feature of our 
studies is a combination of nonperturbative methods applied to a general form of the 
theory, namely, the QCD coupled to QED. Our conclusion is that there can exsist a 
light (nearly massless) axion-like component of the quark-gluon plasma for the given 
range of temperatures. This component accommodates all the chiral and anomalous 
properties of the plasma and is responsible for a plenty of hypothetical effects leading 
to the local V- and CP-violation in the strong interaction. The rest of the matter 
content of QGP obeys the hydrodynamic equations of a nearly ideal fluid. Both of 
the components together form some kind of a superfluid, which we call the "chiral 
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superfluid". This term should be understood in a relative sense as a two-component 
fluid with independent curl-free and normal motions. The separation between two mo- 
tions is provided by a gap in the Dirac spectrum observed in lattice QCD simulations 
with massless quarks. It is worth to mention, that in our case we do not obtain the 
light "superfluid" component as a Goldstone field of a broken continuous symmetry. 
Instead, it appears as a consequence of a nontrivial underlying vacuum structure. At 
the same time, even not being a conventional superfluid, our system reproduces some 
usual properties, such as a (pseudo-) scalar nature of the "condensate", zero entropy 
of the "condensate", Josephson-type relation, etc. Taking into account only regular 
configurations of the "condensate" we reserve the vortex-like solutions for the further 
studies. As a crucial test of our model we propose some phenomenological effects (Sec- 
tion 3.2), which can be proven experimentally, i.e. a response of QGP to the presence of 
strong electromagnetic fields. A concrete description of the experimental consequences 
is currently in progress. 
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